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In the present work, we attempt a brief summary of various areas where nonlinear waves have been 
emerging in the phenomenology of lattice dynamical systems. These areas include nonlinear optics, 
atomic physics, mechanical systems, electrical lattices, nonlinear metamaterials, plasma dynamics and 
granular crystals. We give some of the recent developments in each one of these areas and speculate on 
some of the potentially interesting directions for future study. 

Keywords: 
1. Introduction 

Over the last two decades, there has been an explosion of interest on the dynamics of nonlinear waves 
in lattices. It can be argued that thi s field was actually initiated by the seminal investigation of Fermi, 
Pasta and Ulam (FPU) presented in iFermi et al.\ (1 1954b . who posed the following question: How long 
does it take for long-wavelength oscillations to transfer their energy into an equilibrium distribution in a 
one-dimensional string of nonlinearly interacting particles? This question spurred the activity of a wide 
array of fields including soliton theory, discrete lattice dynamics and KAM theory, all of which remain 
active research fields today. 

The de velopment of soliton t heory , while focusing at the continuum level at first through the inves- 
tigations of IZabuskv & Kruskall d 19671) and the concomitant developments of the theory of integrability 
for the Korteweg-de Vries equation, also offered some of the first particularly interesting examples of 
lattice nonlinear dy namical systems. Among the pri me examples thereof, one can classify the Toda lat- 
tice of Toda (1981), the Ablowitz-Ladi k equation of Ablowitz & Ladik| ( 19761) and the Calogero-Moser 
Af-body problem, see e.g. the works of Calogero ( 197 lb and Moserh 1975 ). 

It was at around the same time that some of the early fundamental suggestions of the role of dis- 
crete solitons in physica l systems were m ade in the 1970s. Some of the pioneers of these contribu tions 
were Davydoy, see e.g. iDavydovl (119721) and Heeger and collaborators, see e.g. ISu et al. I dl979h and 
Heeger et al] ([1988). Davydov proposed the discrete soliton as a tool for understanding energy transfer 
in proteins and their a-helices, while Heeger and collaborators put forth solitonic models for understand- 
ing neutral and charge transport in conducting polymer chains such as polyacetylene and polythiophene. 

However, a true explosion of the field came about in the late 1980s when not only majo r devel - 



opments arose in th e context of the so - called disc rete self-trapping eq uation of lEilbeck et al.\ (119851) 
but also the work of Sievers & Takenol dl988l) and lTakenoeTafl dl988b brought forth the theme of in- 



trinsic localized modes (or discrete breather s) in anharmonic lattic e s. Al thou gh this theme h ad been 
explored earlier e . g.. by lOvchinnikov dl97d). iKosevich & Kovalevl (11974b and iDolgovl (11986b . it was 
Sievers & Takenol ( 1988 ) and Takenoef al. dl988 ) that motivated a number of other researchers to ex- 
plore the existence and stability of such modes more thoroughly and to start obtaining fundamental 
results about their properties, as well as to devise special limits (such as the anti-continuum (AC) limit 
of lattices with uncoupled sites) which could be used to showcase the generic nature of these modes. 
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A landmark in connection to these efforts was the work of iMacKav & Aubrvl ( 1 19941) which established 
rigorously the existence/robustness of such localized modes starting from the AC limit, under minimal 
appropriate (non-resonance) assumptions for chains of nonlinear oscillators. 

Since these early steps, there has been a tremendous amount of developments in the area of nonlinear 
waves in dynamical lattices. This activity has chiefly been fueled by the experimental observation of 
these modes in a wide range of physical systems in biophysics, solid state physics, nonlinear optics, 
atomic physics, granular crystals, and plasma physics, in addition to numerous developments in the 
more classical fields of mechanical and electrical lattices [and even the latter have seen interesting 
recent developments e.g. with the emergence of the field of nonlinear metamaterials, among others, as 
discussed in section |4~2l below 1 . A (partial) list of key physical systems where explanations in terms of 
nonlinear waves in lattices have had an impact includes: 

• The observation of discrete breat hers in complex electr onic materials such as halide-bridged tran- 
sition metal complexes as e.g. in lSwanson et al. dl999h . 



The formation of denaturation bubbles in the DNA double strand dynamics summarized e.g. in 



Pevrardl d2004). 



The emergence of dynamical instabilities, discrete quantum self-trapping and localized modes in 
ul tracold Bose-Einstein conde nsate s in the presence of deep op tical lattices, see e.g. the reviews 
of lBrazhnvi & Konotopl d2004 and lMorsch & Oberthalei1d2006l) . 



The illustration of localized modes of soli tonic, vortex, ring, multipole, surface, gap, necklace and 
numerous other types in the nonlinear optics of evanescently coupled waveguide arrays, as well 



as in th at of biased photorefractive crystals; a relevant recent review can be found in lLederer et al. 
d2008l) . 



The prototypical mechanical realizatio n of solitary waves and localized modes in driven microme- 
chanical cantilever arrays as shown in Sato et al. as well as in coupled torsion pendula 

where a recent example is given bv lCuevas et q/.l (12009b . 



• The examination of interesting propert ies and mobility of l ocali zed modes in nonline ar electrical 
lattices and transmission lines; see e.g. Remoissenetl ( 19991) and English et al. ( 2008 ). 



• The presence of nearly compact solita ry waves in granu lar c rystals consistin g of beads with 
Hertzian elastic interactions reviewed in Nesterenko OOP lb and lSen et all d2008l) . 



The creation of discrete breather type structures in layered antiferromagnetic samples such as 
those of a (C^HsNHs^CuCU crystal, upon imposition of a dynamically unstable continuous wave 
(which f alls into modulationa l inst ability and results in the localized waveforms); see the relevant 
works of English et al.\ (I200lb and lEnglish et al.\ (I2003h . 



These experimental developments have, in turn, motivated a huge range of theoretical advances 
in connection to the properties and interactions of these lattice nonlinear waves and on how these are 
modified in comparison to their continuum siblings. In particular, some of the areas that have seen 
intense theoretical activity concern: 



The existence and stability of the localized modes. 



• The dynamics and interactions of multiple such states in nonlinear lattices. 
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• The effect of external potentials (either local "impurity-like" potentials, or more broad externally 
imposed, field potentials) 

• The effect of short-range (local) versus long-range (nonlocal) interactions. 

• The presence of excitation thresholds for such nonlinear states (especially in higher dimensional 
settings). 

• The statistical mechanics and long-time asymptotics of these lattice systems and the role of the 
localized modes in them. 

• The effects of disorder and the role of the interplay between the linear phenomenon of Anderson 
localization and of the nonlinearity-induced type of localization. 

• The interaction with and scattering of phonon-like modes from discrete solitary waves. 

• The modifications imparted on the waveforms of such localized states and the new possibilities 
for ones such (e.g. with embedded vorticity) afforded by higher dimensional settings. 

Many of these topics have been summarized in a number of reviews b y Braun & Kivshar Aubry 



1997), 



(2003), 



Flach & Willisldl998h.lHennig & Tsironisldl998l) . lKevrekidis e7aZI ( l200lh , [Eilbeck & Johansson! 
Flach&Gorbachld2008h" an d more recently in a number of books e.g. about Klein- Gordon non- 



linear l attices by Braun & Kivsharl d2006l) . as well as about nonlinear Schrodinger ones by lKevrekid"isl 
d2009ab . 

As it can be inferred from the shear volume of review articles on the topic (as well as that of full- 
blown books), it is impossible to give an exhaustive description of any particular aspect or application of 
the subject in anything less than a lengthy review or a book. Our aim herein will thus, by necessity, be 
more modest, but in some sense, also more "forward looking". We aim to give a selective presentation 
of some of recent developments of a few among these areas, which in the present author's evaluation 
hold significant potential for future developments in the field of localized modes, nonlinear waves and 
discrete breathers. Obviously, this bears a significant weight of subjective perception and a personal- 
ized viewpoint about which we caution the reader in advance. On the other hand, our hope is to offer a 
number of interesting directions for future study in these areas which could significantly enhance our un- 
derstanding of nonlinear waves in these settings, promoting their state-of-the-art and offering significant 
connections with present or near-future experimental settings and potentially also with applications. 

We will categorize our presentation based on the different models and the corresponding areas of 
interest. Hence, upon introducing our basic notation, the third section will concern itself with localized 
modes in models of the discrete nonlinear Schrodinger type and their applications in nonlinear optics 
and atomic physics. The fourth section will consider models of the Klein-Gordon type, as well as their 
mechanical and electrical (as well as nonlinear metamaterial) applications. The fifth section will focus 
on models of the FPU type arising in the examination of granular crystals with Hertzian (or modified 
Hertzian) interactions, as well as in the case of study of dusty plasmas. Lastly, in section six, we will 
summarize some of our discussion and offer some concluding thoughts. 



2. Notation 

The prototypical dynamical lattices that will be considered in this review will be the discrete nonlinear 
Schrodinger equation (DNLS) in the form: 

iu» = -e(u n+ i + u„-i) - \u„\ 2 u n , (2.1) 
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the nonlinear Klein-Gordon (KG) lattice of the form: 

ti„ = £(m„ + i+w„_i -2u n )-V'(u n ), (2.2) 
and the FPU-type dynamical lattice of the form 

ti n = V'(u n +\ - U„) - V'(u„ - U n -\), (2.3) 

where u„ is the dynamical variable of interest measured at site n, e is the coupling parameter and V a 
potential (that could be onsite for KG lattices or inter-site for FPU ones). 

It should be noted that in what follows, we will not restrict ourselves to the one-dimensional frame- 
work of Eqs. ( 12.1l )-( 12.3b . but motivated by physical applications, we will often consider higher di- 
mensional variants of the models. In their prototypical form, two-dimensional generalizations of Eqs. 
dPTl-lP) read 

itin,m = £( w H+l,m ~t~ l.m ~t~ w n,m+l ~t~ ^n,m— l) \Un,m\ U n m: (2.4) 

for the DNLS lattice, 

tin,m £ (^n+ 1 ,tn ~t~ ^n— 1 ,m ~t~ ^n,m+ 1 ~^~^n,m—l 4w tt m ) V (w H m ), (2.5) 
for the KG case, while for the FPU-type lattices we have 

ti n ,m = V'( u n+l,m ~ u n,m) ~ V \u n .m ~ u n-l,m) + V (u n .m+l ~ u n,m) ~ V {u n ^ m — U n .m-l)- (2.6) 

This, in turn, clearly suggests how to expand the model beyond two dimensions and into the three- or 
arbitrary-dimensional-case, by sequentially adding more lattice directions and linear (in cases such as 
(I2.4b - (l2.51 l) or nonlinear (as in case ( 12.6b ) interactions along them. 

In the first one of these, the field will be complex and we will most often seek stationary (standing 
wave) type solutions of the form u n = exp(//if )v„; is often referred to as the propagation constant in 
nonlinear optics, or the chemical potential in Bose-Einstein condensates and characterizes the frequency 
of such standing wave solutions. These will satisfy the steady state equation: 

(H - \v n \ 2 )v n = e(v„+i + v„_i). (2.7) 

Additionally, one of the funda mental ideas that we will u se in exploring this class of systems will be that 



of the anti-continuum limit of Ma cKav & Aubrvl (119941) . This is the limit where the sites are uncoupled 



with e = 0. In this case, Eq. ( 12.7b is completely solvable v„ = O,± A /7Texp(/0„), where 9„ is a free 
phase parameter for each site. The fundamental question of interest will involve the persistence of these 
solutions for e ^ 0, as well as their stability. The latter will be often explored by means of linear stability 
analysis, upon imposition of a perturbation of the form: 

u n (t) = exp(ijxt) [v„ + 5 (p„exp(Xt) + q* n exp(X*t))} . (2.8) 

[Notice that the * will be generally used in what follows to denote complex conjugation.] Consideration 
of the resulting equation upon substitution of ( 12.8b to d2.ll ). at 0(8) yields the linearization (infinite- 
dimensional matrix) eigenvalue problem for the eigenvalues X and eigenvectors (p n ,q n ). 

In the case of the Klein-Gordon models, we will focus on discrete breather time-periodic solutions 
for different potentials which are either polynomials in the form V(u„) = Y^=2 a i u n or other functions 
such as V(u n ) = 1 — cos(m„) (in the sine-Gordon case). Notice also that we will often use the operator 
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notation 4im„ = u u +\ + m„_i — 2u„ for the discrete Laplacian. For the FPU lattices, we will primarily 

5 11 

concern ourselves with power law potentials and particularly the Hertzian elastic case of V (r„) ~ r„ , 
but we will also briefly refer to other potentials such as the Yukawa one V(r„) ~ exp(— ar„)/r n (for 
r„ > 0) in the context of dusty plasmas. Notice that in the FPU-type chains, we will also refer to the 
so-called strain formulation which introduces the modified field variable r„ = m„_ i — u n and leads to the 
dynamical equation 

r» = V'(r n+1 )+V'(r n ^)-2V'(r n ). (2.9) 



3. DNLS Models in Nonlinear Optics and Atomic Physics 



3 . 1 Recent Developments 

Many of the investigations in these fields stemmed from the relevance of the DNLS model and its vari- 
ants in the nonlinear optics of fabricated AlGaAs waveguide arrays as in lEisenberg et al. as well 
as in that of Bose-Einstein condensat es in sufficiently deep optical la ttices analyzed in Brazhnvi & KonotopI 
d2004 . lMorsch & Oberthaleri d2006l) . and lKevrekidis et q/.l(l2008al) . 

In the former area, a multiplicity of phenomena such as discrete diffraction, Peierls barriers (the 
energetic barrier that a wave needs to overcome to move from one latti ce site to the next), dif frac- 
tion management (the periodic alternation of the diffraction coefficient) in lMorandotti et al. (1999) and 
gap solitons (st r ucture s localized due to nonl inearity in the gap of th e underlying linear spectrum) in 
Mandelik et al.\ (12004b among others [see also iMorandotti et al 1 dl999h l were experimentally observed. 
More recently, fundamental inves tigations also on th e modulational instability arising in such settings 
for simple plane wave solutions in lMeier et al.\ d2004l) . and o n the role of multip le components and pos- 



Meie r et al. d2003l) have been explored, as 



Suntsov et al. (2006) 



sible four-wave-mixing effects arising from their coupling in . 
well as that of interactions of solitary waves with surfaces in I 

A related area where lattice dynamical models, although not directly relevant as the most appro- 
priate description (since the setting is, in principle continuum with a periodic potential), still yield 
accurate predictions regarding the existence and the stability of nonlinear localized modes is that of 
optically induced lattices in photorefractive media such as Stront ium Barium Niobate (SBN). Since 
th e theoretical suggestio n of such lattices in Efremidi s et al. I d2002l). and their experimental realization 
Fleisc her et al.\ (2003), N eshev et al. tool" d iMartin et al.\ (120041) . there has been a tremendous 



m 



growth in the area of nonlinear waves and solitons in such periodic, predominantly two-dimensional, 
lattices. The continuously growing array of structures that have been predicted and experimentally 
obtained in such lattices inclu des (but is not lim ited to) various discrete multi-pulse pattern s, such as 
the discret e dipole sol i tons i n Yang et al. ( 20041). quadrupole solito ns in Yang et al. ( 2004h. necklace 



solitons in Yang et al. I d2005l). soliton stripes in ,Nesh ev et al.\ d2004l). discrete vortices in lNeshev et al 
(2004), and lFleischer et al.\ d2004h . rotary solitons in lWang et al. I d2006l) . surface solitons in lWang et al 



d2007l) . and so on. Such structures have a definite potential to be used as carriers and conduits for data 
transmiss ion and processing, i n the setting of all-optical communication schemes; see a detailed recent 



Lederer et all (2008 ) 



review in 

Lastly, a completely independent and entirely different field where relevant considerations and struc- 
tures, as well as models of the DNLS type arise is that of ultracold atomic gases in the form of Bose- 
Einstein condensates (BECs) under the confinement imposed by so-called optical lattice potentials. 
These are period ic traps produced by counter-propagating laser beams in one, two or even all three di- 
rections, see e.g. iBurger et al] d2005l) . This field has also experienced a treme ndous amount of dev elop- 
ment in the past few years with some of its landmarks being the prediction [in lSmerzi et al.\ d2002l) l and 
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manifestation [inlCataliotti et al. (2003) ] of modulational instabilities, the observation of gap solitons in 
Eiermann et al. of Landau-Zen er tunneling (tunneling between different bands of the periodic 

potential) in iJona-Lasinio et al. (12003b and Bloch osci l lation s (for matter waves subject to combined 
periodic and linear potentials) in Anderson & Kasevichl ( 1998 ) among many other salient features. 

All of the above settings are described at some level of approximation by a continuum nonlinear 
Schrodinger (NLS) equation with a periodic potential. However, when the periodic potential is suffi- 
ciently deep, it was realized that an equivalent description of the s yste m can be given by a non linear 



dynamical lattice as discussed e.g. in lTromb ettoni & Smerz Hd200ll) . and lAbdullaev et all feOOll) . This 



was further elaborated for BEC settings in lAlfimov et al. 00021) [see also lAblowitz et al.\ d2004l) for a 



related discussion in connection to nonlinear optics], whose presentation we briefly follow here. 

The localized wavefunctions at the wells of the periodic potential can be approximated by Wan- 
nier functions, i.e., the Fourier transform of Bloch functions. Given the completeness of the Wan- 
nier basis, any solution of the NLS equation with a periodic potential can be expressed as u(x,t) = 
Ln.a c n , a (t)w n ,a(z)> where n and a label wells and bands of the periodic potential, respectively. 

Substituting the above expres sion into the dynamic al equation, and using the orthonormality of 
the Wannier basis as indicated in lAlfimov et al.\ (120021) . we obtain a set of differential equations for 
the coefficients. Upon suitable decay of the Fourier coefficients and the Wannier functions' prefactors 
(which can be systematically checked for given potential parameters), the model can be reduced to 



.dc na 
l — ; — 
dt 



&0,aC n ,a + &l,a {c n -l,a + 



c n +l,a)+ 

cn,a 2 ,a 3 



•njramn * 



/7,0!3 1 



(3.1) 



1 jjrimimn^ poo 

where W aa \d 2 a 3 = J_. 



w nta w„ liai w„ 2iaz w n3ta} dx is the Wannier function overlap integral and 



k/L 

<»«,a = z— / E a (k)e 



—iknL 



2k J-k/l 



dk , 



(3.2) 



is the Fourier transform of the energy E a (k) of the Bloch mode corresponding to wavenumber k for a 
potential of spatial period L. This degenerates into the so-called tight-binding model 



.dc na 
l — 

dt 



Oh.aCn.a + Sh,a (c n -l,a + C B +l,a) + W 



1111 l L «;a| L n,a, 



(3.3) 



i.e., the DNLS equation, when restricting consideration to the first band. Higher-dimensional versions 
of the latter are of course physically relevant models and have, therefore, been used i n various studies 



concern ing quasi-2D and 3D BECs confined in strong optical l attices as illustrated in iKevrekidis et al 
(2008a) and similarly in 2D waveguide arrays as summarized in IKevrekid is (2009a). 



The relevance of the DNLS model in this wide array of application areas prompted an intense activity 
in investigating the fundamental solitary wave solutions of the model. We now focus on analyzing these 
prototypical solutions and their stability in the context of Eq. ( 12. 11 1. In particular, for stationary solutions, 
upon multiplying Eq. ( 12.7b by v* and subtracting the complex conjugate of the resulting equation leads 
to: 

v*„v„+i - v„v* +1 = const. ^> 2arg(v„ + i ) = 2arg(v„), (3.4) 

for localized solutions (i.e., vanishing as n — > ±°°). Using the scaling freedom of the equation to set jX = 
1 allows us to infer that the only states that will persist for finite £ will be t he one s containing sequences 
with combinations of v„ = ±1 and v„ = 0. The work of lAlfimov et al.\ (12.004ft provided a systematic 
numerical classification of the most fundamental of the resulting sequences and their bifurcations. Some 
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FIG. 1. (Color online) A sampler of solutions of the discrete nonlinear Schrbdinger equation in 1-, 2- and 3-spatial lattice 
dimensions. The top left shows the always stable fundamental single soliton. The top middle and right show the out-of-phase 
two- and three-site excitations (which can be stable near the anti-continuum limit). The bottom left (4-panel set including the 
amplitude, phase, real and imaginary part contour plots of the solution) shows the prototypical vortex of topological charge 5=1 
(which is also stable near the AC limit). The bottom middle panel contains a rhombic S = 2 generalization thereof with 8 sites 
(which may also be stable), while the bottom right contains the also stable near the AC limit 6-site diamond three-dimensional 
configuration. 

of the prototypical examples of sequences involving one, two and three sites are illustrated in Fig. Q] 
While in our discussion below, we restrict our consideration to the case of focusing nonlinearity, our 
results can be extended to the defocusing case [with the opposite sign of the nonlinearity in Eq. (12. 1 I II 
via the so-called staggering transformation w n = (— l)"u n . This converts the defocusing nonlinearity 
into a focusing one, with an appropriate frequency rescaling which can be trivially absorbed in a phase 
or gauge transformation. 

We now turn to the examination of stability. Using the linearization ansatz of Eq. ( |2.8t and a simple 
eigenvector rotation leads to the equivalent symplectic formalism of the linear stability problem /Jzfw = 
Xw, where _£f is a diagonal, self-adjoint matrix with the operators Jzf + and Jr?_ given below in its 
diagonal entries and J is the symplectic matrix. Jzf+ and Jzf- are defined by: 

(1 -3v^)a„-e(fl„ + i+fl„_i) = Jz?+a„ = -Xb n (3.5) 
(l-vl)b n -e{b n+1 + b n -i) = &-b n = Xa n . (3.6) 

Using once again the AC limit, we assume a v n with N "excited" (i.e., ^ 0) sites; it can then be straight- 
forwardly inferred that for e = these sites correspond to eigenvalues X + = —2 for Jz? + and to ones 
with A_ = for Jzf_. In turn, these imply the existence of N eigenvalue pairs with X 2 = for the full 
problem. These N vanishing eigenvalue pairs are potential sources of instability, since N — I of those 
will become nonzero, upon departure from the AC limit, given that there is a single symmetry left for 
nonzero e, namely the U(l) invariance i.e., the invariance with respect to phase. The key issue for 
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stability purposes is to identify the location of these N — I small eigenvalue pairs. One can manipulate 
Eqs. ( 13.5b — (f3~6l > into the form: 



&_b n = -X 2 &7 l b n 



In the vicinity of the AC limit, the effect of «5f + is a multiplicative one (by —2). Hence: 



1 



lim(b„,^+ l b„) = -- => X z = 2y=Ub„,#-bn) 
£->0 2 



(3.7) 



(3.8) 



In light of this calculation, the full problem becomes equivalent to the determination of the spectrum 
of Jz?_ . A crucial fact in that regard is that v„ is an eigenfunction of Jzf_ with A_ = 0. Then, direct 



use of the Sturm comparison theorem for difference operators [see e.g. lLevv & Lessmanl dl992l) l leads 
us to conclude that if the number of sign changes in the solution at the AC limit is m (i.e., the number 
of times that adjacent to a +1 is a —1 and next to a —1 is a +1), then n(J?f_) = m and therefore 
from Eq. ( 13.8b . the number of imaginary eigenvalues pairs of J££ is m. This, in turn, implies that the 
number of real eigenvalue pairs is consequently (N — 1 ) — m. Another immediate conclusion is that 
unless m =N—l, i.e., unless adjacent sites are out-of-phase with respect to each other, the solution 
will be immediately uns table upon departure from the AC limit. This analysis originally presented in 
Pelinovskv et al.\ (12005 ah is also consistent with the general eigenvalue count of iKapitula et al.\ d2005l) . 
It should be noted that these m imaginary eigenv alue pairs end up possessing, so-called, negative Krein 
signature as analyzed in lPelinovsky et al. I d2005ah . This implies that upon collision with other eigenvalue 
pairs, such as the ones stemming from the continuous spectrum, they will lead to instability through a so- 
called Hamiltonian-Hopf bifurcation and a complex eigenvalue quartet. In this problem, the continuous 
spectrum eigenmodes of eigenfrequency of wavenumber k are characterized by the dispersion relation 
A = ±i(l +4£sin 2 (£/2)), and hence correspond to a band of width 4e, separated by a unit distance 
from the origin along the imaginary axis of the spectral plane (X r , A;) of the eigenvalues X = X, - + iX;. 

Equation (13. 8t can also be used in a quantitative fashion to identify the relevant eigenvalues pertur- 
batively for the full problem. In particular, this amounts to considering the eigenvalues of _§f_ emanating 
from through the perturbed eigenvalue problem: 



J? W bl?\ 



(3.9) 



where Jzf_ = Jz?l 0) + e^i' + @{e 2 ) and a similar expansion has been used for the eigenvector b n . For 
these eigenvalues, we know that A_ = eyi + ^(e 2 ). Projecting the above equation to all the eigenvectors 
of zero eigenvalue of _£f_ , one can explicitly convert Eq . d3.91 > into an eigenvalue problem of the 



form Mc = J\c, as was shown in IPelinovskv et cd. ( 2005al) . The matrix M has off-diagonal entries: 
M n:n +i =M n+ \ „ = — cos(0„ + i — 0„) and diagonal entries M„„ = cos(0„_i — 6„) +cos(0„ + i — 6„). The 
subsequent computation of the leading order correction yi (t° the eigenvalues of ££J) allows us to 
calculate the perturbed eigenvalues of the full problem A = i^/Zeyf. 

Let us consider as a case example the one-dimensional configurations with two-adjacent sites with 
phases Q\ and Ot. Then, the matrix M becomes: 



M = 



cos(0! -0 2 ) 

-cos(ei-e 2 ) 



-cos(0i-0 2 ) 
cos(0! - 2 ) 



(3.10) 



which leads to A 2 = and A 2 = 4ecos(0i — 2 ). Notice that, as expected, for same phase excitations 
(01 = 2 ), the configuration is unstable due to a real eigenvalue pair, while the opposite is true if Q\ = 
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2 ± 7t. This type of analysis is possible e.g. for 3-site configurations with phases d\ 23 (or for that 
matter, for arbitrary numbers of excited sites). In the 3-site case, one of the eigenvalues of M is again 
(this is true for any configuration due to the U(l ) invariance), while the other two are given by: 



7i 



cos(0 2 -0i)+cos(0 3 - 2 ) 



± WCOS 2 (0 2 - 01 ) - COS (0 2 - 01 ) COS (0 3 - 2 ) + COS 2 (0 3 - 02 )■ 



(3.11) 



This formulation of the existence and stability pr oble ms has been generalized to different settings, 
such as higher dimensio ns in Pelinovsky et al. (2005b), and lPelinovskv et al.\ ( 120081) or multi-component 
systems analyzed e.g. in lKevrekidis & Pelinovsky (E006l) . Arguably, the principal difference that arises 
in the higher dimensional settings is that the wave profile may contain sites excited over a contour (as 
opposed to along a straight line). In that case, for the N excited sites around the contour, the persistence 
(Lyapunov-Schmidt) con diti ons can be obtained as a generalization of Eq. d2.lt that reads [see e.g. 
Pelinovsky et al\ d2005bb and lAlexander et al.\ d2004l) l: 



sin(0i-0 2 ) = sin(02 - 3 ) 



sin(0jv — 0i ). 



(3.12) 



This highlights a key difference of these higher dimensional settings, namely that not only "solitary 
wave" structures with p hases € |0,7r| are poss ible, but also both symmet ric and asymmetric vor- 
tex families presented in iPelinovskv et al. d2005bl) and lAlexander et al.\ ( 120041) may, in principle, exist 
[a note of caution, however, is that Eq. ( 13.121 ) provides only the leading order persistence condition 
and one would need to also verify t he corresponding conditio ns to higher order to confirm that such 
solutions persist as was dis cussed in Pelinovsky et al. ( 2005bl)l. These vortex solutions had been pre- 
dicted numerically earlier in lMalomed & Kevrekidisl(12001l) and lKevrekidis et al. Id2004al) and have been 
observed exper i menta lly in the optical setting of photorefractive crystals by iNeshev et al] (120041) and 
Fleischer et al.\ (120041) . The stability analysis is also possible for these higher dimensional structures, 
although the relevant calculations are technically considerably more involved. In fact, these types of 
calculations are possible not only for square lattices, as in the cases highlighted above, bu t for more 



compl ex lattice settings as well including hexagonal and honeycomb lattices as shown in lLaw et al. 



(2008). A number of interesting conclusions arise therein including e.g. the instability of a lower topo- 
logical charge (with vorticity 5=1) vortex and the stability of a higher topological charge on e (with 
S = 2); such results have been recently also confirmed experimentally by iTerhalle et al. ( 2009 ). Some 
case examples of interesting two-dimensional structures are given in Fig. Q] including a prototypical 
vortex cross of topological charge S = 1, as well as a potentially stable S = 2 generalization t hereof in a 
rhombic configuration discussed recently in lOster & Johansson d2006l) and lLaw et all d2009l) . 

However, it should also be noted that the relevant theory can be formulated in an entirely general 
manner: we give the outline below, as well as illustrate some prototypical higher-dimensional (i.e., 3D) 
structures. In the multi-dimensional case, the problem of existence of stationary standing wave solutions 
can be formulated through the vanishing of the vector field F„ of the form: 



f„(<M) 



(i-|</>„| 2 )4>„-£i:4>„ 

(l-\<j> n \ 2 )<j>*-eE<j>* 



(3.13) 
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Upon defining the matrix operator: 

/ 1-2|(^„| 2 



¥P — 




1 
1 



(3.14) 



where the s± e . denote the shift operators along the respective directions and e, are the corresponding 
unit vectors in each of the three lattice directions, the stability problem reads affy = ;A l//. Here each 
2x2 bloc k of the matrix a is the diagonal matrix with elements (1,-1) along the diagonal for each 
node «; see lPelinovskv et al.\ (120081) for a detailed discussion. The existence problem, on the other hand, 
is connected to 3tf through: Jf? = DaF(0,O). At the AC limit of e = 



(^ (0) )„ 



1 
1 



neS J 



-1 

,-2i6 



-1 



n E S, 



where S is the set of excited sites. Then the eigenvectors of zero eigenvalue will be of the form: 



It is helpful to define the projection operator: 



(f2)„), n£S, 



and to decompose the solution as: 



<j) = <j) (0) (G) + (peX. 



(3.15) 



(3.16) 



Then, one can obtain the Lyapunov-Schmidt persistence conditions analyzed in lPelinovsky et al.\ 12008) 
as: 

g(0,e) = ^F(> (o) (0) + <p(0,£),£) =0. (3.17) 

As proved in iPelinovskv et al. | d2008l) . the configuration 0^(9) can then be continued to the domain 
£ € G(0) (i.e., for couplings in the neighborhood of the AC limit) if and only if there exists a root 0* 
of the vector field g(0,e). Moreover, if the root 0* is analytic in e G ^(0) and 0* = 9q + &{e), the 
solution of the difference equation is analytic in e £ ^(0), such that 



= </> (0) (9*) + <p(0*,£) = (O) (0 O ) + E e^ w (0 o ). 

k=\ 



(3.18) 



Within the same formulation, one can establish a general stability theory, provided that the relevant 
solution persists for £ ^ 0. If the operator 34? has a small eigenvalue of multiplicity d, such that 



then the full Hamiltonian eigenvalue problem admits (2d) small eigenvalues A. 



These are such that A = £ k ^ 2 ^k/2 + ^(e^ 2+1 ), where non-zero values X k / 2 are found from 



*/2 



odd A:: ^ {k) a = ^Xl /2 a, 



even 



k: JZ {k) a + \\ k /2^ {k) a = X - tf /2 a 



(3.19) 



(3.20) 
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where 



= D 0g «(0 o ), 



and k' = {k— l)/2. For more details, we refer the interested reader to IPelinovsky et al.\ (12008b . A 



prototypical example of a configuration (the so-called discrete diamond) that the above theory predicts to 
exist and be dynamically robust near the AC li mit in 3D is illustrated in the bottom rig ht panel of Fig, [p 
Other s uch configurations can be found e.g. in IPelinovsky et al.\ 020081) . as well as in iKevrekidis et al. 
d2004bl) andlCarretero-Gonzalez et all d2005l). 



3.2 Future Perspectives 

Although, based on these recent results, much is presently understood about the statics and stability of 
DNLS type systems, a number of themes still remain relatively unexplored. We now present a few of 
the themes that we believe hold significant interest for future studies. 



3.2.1 A general theory and applications of long-range interactions. While the setting of local (i.e., 
nearest-neighbor) interactions has been substantially explored in the DNLS case, this can far less be 
argued to be the case for the setting of non-nearest neighbor effects. Some early investigations of 
one dimensional settings identified the role of non-nearest neighbor interactions in creating bistability 
of the fundamental soliton solutions and leadi ng to interesting possibi lities for switching among sta - 
ble branches (upon suitable perturbations) as in lJohansson et al.\ ( 11998b and lChristiansen et al] (11998b . 
However, a more systematic investigation of existence and stability issues, especially for more compli- 
cated solutions was not offered. 

More recently, a number of physically-minded works have argued the relevance of the inclusion 
of such non-nearest neighbor interactions i n a variety of settings, involving p redominantly optical 
waveguide arrays. For instance, the work of lEfremidis & Christodoulidesl (12002b argued that a zigzag 
waveguide array could be considered as a quasi-one-dimensional chain in which the relative strength 
of nearest-neighbor and next-nearest-neighbor interactions can be tuned on the basis of the relevant 
geometry. Corresponding ideas for involving non-nearest-neighbor interactions have also been gener- 
alized to two-dimen sional settings. Although limited to the linear propagation at least in the realm of 
Szameit et al.\ (12008b . even in that setting, there is interesting phenomenology arising from the competi- 
tion of the different types of interactions, leading potentially to diffraction-free pro pagation for suitable 
wavenumbers in the center of the Brillouin zone. On the nonlinear side, the work of lKevrekidisI d2009bl) 
indicates that both the existence problem and the stability may be crucially modified with respect to 
the standard DNLS case by the presence of nonlocality. A typical example of the former involves the 
emergence of ID solutions that have arbitrary (i.e., different than 0, %) phases. A typical example of the 
latter involves t he stabilization of unstable structures [such as a prototypical class of S = 2 vortices in a 
square lattice of IPelinovsky et al. d2005bb l. It should also be noted that another area that may enhance 
the relevance of long-range interaction considerations is that of dipolar Bose-Einstein con densates (such 



as 52 C r) in the presence of optical lattices in atomic physics; see e.g. the recent review of lLahaye et al 



d2009b 



These directions and findings seem to warrant a more systematic investigation of the effects of 
nonlocality, ideally as a general function of properties (e.g. the "interaction range") of the kernel. Such 
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a general setting could be: 

N 

iu„ = -e £ a nm u m - \u„\ 2 ii„ (3.21) 

m=l 

for different types of interaction kernels a nm (although the nonlocality could also in principle, or addi- 
tionally be imposed on the nonlinear term). It would appear to be timely and relevant to explore how 
existence, stability and dynamics of solitary waves are progressively modified, as we depart from the 
well-understood local limit. 



3.2.2 Intersite Lattices, Discretizations, Symmetries and Traveling Waves. One of the directions that 
also hold promise in the setting of DNLS models is that of different types of discretizations and their 
properties in comparison to the more fundamentally well-understood model with the centered-difference 
Laplacian and the onsit e cubic nonlinear i ty. As a motivating example in this direction, we briefly dis- 
cuss the recent report o flKevrekidis et al\ (I2009cl) which examined the two-dimensional Ablowitz-Ladik 
model of the form: 



£ 1 ,m Un— 1 ,m ~t~ Mn,m+ 1 ~t~ w /7,m— 1 ^Un.ni) 

— \u\ 2 (u„ + \ Jn + U„-\ m + ll„j n +\ + U„,m-l) ■ 



(3.22) 



The fundamental difference of this, so-called AL-NLS, model from the standard DNLS is that in ad- 
dition to the centered-difference approximation of the Laplacian, a nearest-neighbor average is used to 
discretize the cubic nonlinearity of the continuum limit a\u\ 2 u (instead of a local term C7\ii„. m \ 2 u„ M in 
the DNLS). In the one-dimensional se tting, such a centered-diff erence discretization of the nonlinearity 
is completely integrable, as shown bv lAblowitz & Ladik il976l) . giving rise to exact hyperbolic-secant 
solitonic solutions that can travel at arbitrary speeds. This well-known (since the 1970s) fact already 
showcases the special properties that can emerge upon dif ferent types of discretization. 



However, the recent work of Kevrekidis et al. 



(I2009d) illustrated that surprising features may arise 
because of such discretizations in higher-dimensional settings as well. In particular, in the 2D case, 
it was found that the DNLS solitons become unstable at some critical threshold (e.g. of the coupling 
strength e) as the continuum limit is approached (e.g. by increasing e). The solitons then remain unsta- 
ble all the way to the continuum limit (of e — > °°) with the relevant instability eigenvalue approaching the 
origin of the spectral plane. This is because the 2D continuum NLS model is critical i.e., it is marginally 
unstable with respect to collapse . In the critical case, the continuum model is invariant under rescaling, 
see e.g. the relevant analysis of ISulem & Sulem dl999h . On the other hand, the AL-NLS approaches 



this continuum limit in a completely different way; in particular, while the AL-NLS two-dimensional 
solitons become unstable within a narrow range of parameters (such as e), they subsequently become 
restabilized and remain stable as the continuum limit is approached. This suggests the remarkable fact 
that while the two models (DNLS and AL-NLS) are identical at the limit, infinitesimally close to the 
limit, their dynamics is substantially different. While the DNLS solitons are exponentially (yet weakly) 
unstable, the AL-NLS ones are dynamically stable. 

Further consideration of this feature indicates that it is a particular trait of critical settings, which 
are at the very special separatrix between subcritical settings where the solitary waves are dynamically 
stable and supercritical ones, where the waves are exponentially unstable. In this critical case, the linear 
spectrum possesses an additional zero eigenvalue pair (associated with the pseudo-conformal invari- 
ance), which permits the reshaping of the solution under the action of the group of rescalings, and hence 
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paves the way for the emergence of self-similar collapse. Discreteness can then shift this pair along the 
imaginary axis or along the real axis. The AL-NLS discretization turns out to be a prototypical example 
whereby the eigenvalue pair formerly associated with the pseudo-conformal invariance is perturbed in 
a stable way (moves along the imaginary axis of the spectral plane), upon discretization and hence this 
model allows infmitesimally small spacings to give rise to collapse-free dynamics. 

The above discussion of issues pertaining to symmetry raises the more general question of devis- 
ing different types of discretizations and examining their symmetry properties and the connection of 
these to the dynamical features of solitary waves. Assuming that discretizations of the continuum NLS 
model will respect the phase/gauge-invariance of that model, another key symmetry whose impact on 
discretizations has been examined fair ly extensively recently is the invariance with respec t to transla- 
tions Motivated by the early work of ISpeight & Ward! d 19941) . ISpeightl d 19991) and that of iKevrekidisI 
( 12003b . there has been a considerable volume of literature developed recently on the subject of suitable 
discretizations of Klein-Gordon, as well as of NLS m odels which preser ve (for their stationary solu- 
tions) the property of translational invariance; see e.g.. |Pelinovskyl (120060 and references therein. It is 
well-known that such an invariance is absent in the standard DNLS model which admits single-humped 
discrete solitary waves which can only be centered on a site or half-way through between two sites (but 
not other such types of waveforms, contrary to what is the case with the translationally invariant models 
where the solitar y waves can be cen tered anywhere). In this context, and for the NLS-model discretiza- 
tions the work of P elinovskv d2006l) unified some of the earlier findings and methods, illustrating that 
the most general translationally invariant discretization of the cubic model is of the form: 



-e(w„+i +m„_i -2u n ) - 2r\)\u„\ (u n+ i +m„-i) - %u n (u* n+x +u„-i) 

-<•; [( 2 I««I 2 + K+il 2 + K-i l 2 )"» + + ««+i^-iH + (« 2 +i + m 2 -iK] 

-7?(K+i | 2 + K-i | 2 )(m„+i + m„-i) - V [ul +i u*_ l + ul^u* +l ■ 



(3.23) 



\Un+\\ 2 Un-\ — |M)i-l| 2 Mn+ll , 



where (x, ^ , rj, v) are real-valued parameters. When all four parameters vanish, one recovers the AL- 
NLS limit. 

This discussion, in turn, leads to the interesting additional question of the possibility of identifying 
traveling wave solutions in DNLS (and related) lattices. This question was considerably controversial 
during the previous decade with many contradictory results and conjec tures; see e.g. the discu s sion i n 
Melvin et al. ( 2006h and references therein. Starting from the work of Gomez-Gardenes et al.\ d2004l) . 
which illustrated numerically that in DNLS such traveling waves cannot be localized, this question 
was subsequently answered bot h in the DNLS setting a nd in that of more complex models such as the 
saturable nonlinear proposed in lHadzievski et al.\ (120041) in the form: 



iu n 



-£A 2 u„ 



1 



(3.24) 



What the direct numerical simulations of Hadzievski et al. 



roborated numerically in the work of Melvin et al. ( 2006) . iMelvin et al.\ d2008l) and quasi-analytically 



d2004l) : nitially sugge sted and which was cor- 



through the asymptotic expansions of lOxtobv & Barashenkovl (2007) was that indeed in DNLS, the res- 
onance of such traveling wave ("embedded soliton") type solutions with the (modified, for traveling 
solutions) continuous spectrum gave rise to "nanoptera" i.e., solutions with a non-vanishing tail. How- 
ever, in different types of nonlinearities such as the saturable or cubic-quintic ones, these works showed 
that it is possible to make the prefactor of such a resonance (the so-called Stokes constant) vanish, 
thereby producing non-generic, yet exact, exponentially localized traveling solutions in these models. 
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The above findings, in turn, beg the somewhat open-ended question: is there some simple diagnostic 



which ma y indicate the existen ce of such traveling sol utions in a model? The works of Hadzievski et al. 
( 2004 ) and Melvin et al. ( 20061) [see also the studies of Carretero-Gonzalez et al. ( 20061) and Chong et al. 



(2009) for corresponding discussions of the cubic-quintic discrete model] hinted at the use of the so- 
called Peierls-Nabarro barrier, i.e., the energy difference between on-site and inter-site centered so- 
lutions. In particular, they argued that vanishings of this barrier should be physically expected to be 
connected with the possibility of such traveling solutions. However, this connection has not been made 
rigorous, as of yet. More generally, exploring systematically the connection of Peierls-Nabarro barriers, 
the existence of traveling solutions, the calculation of the Stokes constant, and the potential for under- 
lying symmetries (including a possible semi-discrete type of Galilean invariance) could be extremely 
interesting topics for further exploration over the next decade. Understanding such properties either on 
a model-specific basis, or, ideally, based on more general/fundamental principles is an important open 
direction for these nonlinear d ynamical lattices. At this point, we should also mention in passing the 
very interesting recent work of Barashenkov & van Heerdenl ( 2008 ). There, the examination of excep- 
tional discretizatio ns [ bearing an effective transla tional invariance through the map type approach of 
Pelinovskv (2006) and Barashenkov et al. d2005h l led to some ingenious suggestions on how to dis- 
cretize so as to preserve genuinely traveling localized excitations i.e., kinks in the discrete sine-Gordon 
and discrete 4 settings. It would be extremely worthwhile for such methodologies to be brought to bear 
also in DNLS and related settings. 



3.2.3 Statistical Mechanics of the DNLS and Approach to Equilibrium. One of the questions that has 
received a somewhat limited amount of attention is that of the statistical mechanics of the DNLS model 
and how a given initial condition approaches its asymptotic (equilibrium or ne ar-equilibrium) dynamics . 
One of the early studies of such questions for the DNLS lattice appeared in iRasmussen et al. d2000l) . 
There, a regime in phase space was identified wherein regular statistical mechanics considerations apply, 
and hence, thermalization was observed numerically and explored analytically using regular, grand- 
canonical, Gibbsian equilibrium measures. However, the nonlinear dynamics of the problem renders 
permissible the realization of regimes of phase space which would formally correspond to "negative 
temperatures" in the sense of statistical mechanics. The novel feature of these states was found to 
be that the energy spontaneously localizes in certain lattice sites forming breather-like excitations (as 
observed numerically and experimentally). Returning to statistical mechanics, such realizations are 
not possible [since the Hamilt onian is unbounded, as is seen by a simple scaling argument similar to 
the continuum case studied in iLebowitz et al.\ (119881) 1 unless the grand-can onical Gibbsian measure 
is refined to correct for the unboundedness. This correction was argued in Rasmussen et a/1 d2000h 



to produce a discontinuity in the partition function signaling a phase transition which was identified 
numerically by the appearance of discrete breathers. 

More recently, the statistical understanding of the formation of localized state s an d of the asymptotic 



dynam ics of the DNLS equation has been examined in the w o rks o f Rumpf ( 2004)) and lJohansson & Rasmussen 
(120041) . Using small-amplitude initial conditions, iRumpfi (120041) argued that the phase space of the 



system can be divided roughly into two weakly interacting domains, one corresponding to the low- 
amplitude fluctuations (linear or phonon modes), while the other consists of the large-amplitude, local- 
ized mode nonlinear excitations. Then, based on a simple partition of the energy H = H +H + and 
of the norm P = P +P+, into these two broadly (and also somewhat loosely) defined fractions, one 
smaller than a critical threshold (denoted by "— ") and one larger than a critical threshold (denoted by 
"+")> it is possible to compute thermodynamic qua ntities such as the entropy in this localization regime. 
In particular, one of the key results of the work of lRumpfl (12004b is that, for a partition of K sites with 



4- 
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large amplitude excitations and N — K sites with small amplitude ones, an expression is derived for the 
total entropy (upon computing S-, S+ and a permutation entropy due to the different potential location 
of the K and N — K sites). This expression reads: 



S = N\n£2 



2E 4 



-inr, 



(3.25) 



where Q = (AP 2 - E 2 ) / '(4A_(JV - K)) and T = 2P-N/P+E-, while K = P 2 /(2E + ). While some 
somewhat artificial assumptions are needed to arrive at the result of Eq. (13.25t (such as the existence of 
a cutoff amplitude radius R in phase space), nevertheless, the result provides a transparent physical un- 
derstanding of the localization process. The contributions to the entropy stem from the fluctuations [first 
term in Eq. (I3.25H and from the high amplitude peaks (second term in the equation). However, typically 
the contribution of the latter in the entropy is negligible, while they can absorb high amounts of energy. 
The underlying premise is that the system seeks to maximize its entropy by allocating the ideal amount 
of energy H to the fluctuations. Starting from an initial energy H , this energy is decreased in favor 
of localized peaks (which contribute very little to the entropy). The entropy would then be maximized 
if eventually a single peak was formed, absorbing a very large fraction of the energy while consuming 
very few particles. Nevertheless, practically, this regime is not reached "experimentally" (i.e., in the 
simulations). This is because of the inherent discreteness of the system which leads to a pinning effect 
of large amplitude excitations which cannot move (and, hence, cannot eventuall y merge into a single 
one) within the lattice. Secondly, the growth of the individual peaks, as argued in IRumpfl ( 120041) . stops 
when the entropy gain due to energy transfer to the peaks is ba lanced by the entropy loss due to transfer 
of power. While placing the considerations of IRumpfl ( 120041) in a more rigorous setting is a task that 
remains open for future considerations, this conceptual framework offers considerable potential for un- 
derstanding the (in this case argued to be infinite, rather than negative, temperature) thermal equilibrium 
state of coexisting large-amplitude localized excitations and small-a mplitude background fluctuations. 

On the oth er hand, the work of Johansson & Rasmussen! (|2004) extended the considerations of the 
earlier work of R asmussen et al. d2000l) to the generalized DNLS model of the form: 



iu„ + (u n+ \ +u n -\) + \u„\ 2a u„ = 



(3.26) 



with a being a free parameter within the nonlinearity exponent. Furthermore, in the work o flJohansson & Rasmussen 
(|2004l) . the connection of these DNLS considerations with the generally more complicated Klein- 
Gordo n (KG) models was discussed. Much of the above mentioned phenomenology, as argued in lRumpf 
( 2004), is critically particular to NLS type models, due to the presence of the second conserved quantity, 
namely of th e I 2 norm; this feature is absent in the KG lattices, where typically only the Hamiltonian 
is conserved. IJohansson & Rasmussen! d2004l) formalize the connection of DNLS with the KG lattices, 
by using the approximation of the latter via the former through a Fourier expansion whose coefficients 
satisfy the DNLS up to controllable corrections. Within this approximation, they connect the conserved 
quantity of the KG model to the ones of the DNLS model approximately reconstructing the relevant 
transition (to formation of localized states) criterion discussed above. However, in the KG setting this 
only provides a guideline for the discrete breather formation process, as the conservation of the norm is 
no longer a true but merely an approxima te conservation law. This is observed in the dynamical simu- 
lations of lJohansson & Rasmussenl (120041) . where although as the amplitude remains small throughout 
the lattice the process is well described by the DNLS formulation, when the discrete breathers of the 
KG problem grow, they violate the validity of the DNLS approximation and of the norm conservation; 
thus, a description of the asymptotic state and of the thermodynamics of such lattices requires further 
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elucidation that necessitates a different approach. This is another interesting and important problem for 
future studies. 



It is also worth pointing out that many of t he considerations of works such as that o flRasmussen et al. 



d2000h or of lJohansson & Rasmussenl (120041) are intrinsically one-dimensional in nature and can not be 
straightforwardly generalized (per the nature of the transfer integral technique used) to higher dimen- 
sions. Hence, the issue of the asymptotic dynamics is perhaps more pressing in such higher-dimensional 
settings, especially given their connections to collapse (sufficiently close to the continuum limit). But 
even in one- dimension, numerous que stions remain. One such interesting question was raised by the re- 
cent work of lKevrekidis et a/.l(l2008bl) . This work offers an understanding of the fundamental threshold 
(for initial data supported on a single site) between discrete dispersion for sub-critical initial amplitude 
and nonlinearity-induced localization, for sufficiently high initial amplitude. However, it also poses the 
question: once we know localization will ensue, which "member" of the mono-parametric family of lo- 
calized solutions will the dynamics "select" for the asymptotic relaxation state? Numerical experiments 
illustrate that the initial condition will shed both some mass and some energy, eventually converging to 
one particular member of that family of stationary solutions. But what is missing is the guiding prin- 
ciple of such a selection. Interestingly, while the single-site supported initial data can give localization 
in the DNLS model (when it possesses a sufficiently large amplitude), in the AL-NLS model, such data 
always leads to dispersion. In the latter case, the machinery of the theory of integrability can be brought 
to bear to appreciate the effect of dif ferent types of compactly supported initial conditions as shown in 
Espinola-Rocha & Kevrekidis! (120091) . 



3.2.4 Interplay of Nonlinearity and Disorder and its Implications on Anderson Localization. One of 
the particularly interesting recent developments on the front of DNLS equations has been the examina- 
tion of the interplay of nonlinearity with disorder, in an effort to explore how the presence of the former 
affects the funda mental phenomenon of Anderson localization due to the latter as first illustrated in 
Andersonl (119581) . This topic has spurred a significant controversy within the physics and mathematics 
communities, given its fundamental nature and the somewhat conflicting result s reported recently. In 



part icular, a number of recent experimental studies in opt ical [see the works of | Schwartz et al.\ (J2007) 
and lLahini et al.\ d2008l) l and atomic [see e.g. the work of ISanchez-Palencia et al.\ (12007b 1 physics set- 
tings have reported the observation of the linearly induced exponential localization proposed by Ander- 
son for a disordered linear lattice. Nevertheless, in a series of recent computational publications in the 
p hysics literature by at least two separate grou ps ( whose theoretical arguments produ ce different results) 



Kopidakis et a/.l (120081) . iFlach et al. I (l2009h . and lPikovsky & Shepelyan sky (2008), it has been argued 



m 

that nonlinearity essentially destroys Anderson localization by producing a subdiffusive scaling of a 
quantity such as the second momentM? = £(« — no) 2 \u n (t)\ 2 (where |«„(f)| 2 plays the role of the wave- 
function density in the Schrodinger problems considered therein). This subdiffusive scaling implies that 
progressively more distant lattice sites are occupied (although slowly) and the process is not stopped 
by the trappin g anticipated on the basis of the Anderson mechanism. More specifically, it has recently 
been argued in lFlach et al.\ (12009b that there exist three regimes within the system's spectral dynamics. 
For weak nonlinearity (smaller than necessary to induce transitions between linear modes), the system 
is in a transient Anderson localization regime, but eventually it detraps from it and grows according 
to the subdiffusive scaling M2(t) ~ t a , with a < 1. For intermediate nonlinearities, the subdiffusion 
of Mi is initiated immediately, while large nonlinearity leads to localization of a fraction of the initial 
wavepacket, w hile the "radiative" r emainder is also subject to this subdiffusive expansion. In the more 
recent work of Skokos et al. d2009h . the relevant considerations were extended to more complex initial 
data (not only of single site, but also of single mode or finite size), and also a similar phenomenology 
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was proposed in iKrimer et al. d2009l) to arise in the setting of DNLS equations with linear (so-called 
Stark ladder) potentials. In the latter, a discrete, equidistant spectrum is also known to arise (but now due 
to the linear potential), hence the same subdiffusion type dynamics is expected to ensue in the presence 
of nonlinearity. This controversy is extremely interesting, since it highlights the potential "fragility" of 
the Anderson localization regime, indicating that even very weak nonlinearity is sufficient to eliminate 
the relevant phenomenology. 

However, all the above considerations, at present, stay within the realm of (admittedly, very strongly 
suggestive) numerical experiments. Neither real experiments (e.g., with optical waveguides or BECs in 
optical lattices) or rigorous mathematical arguments have been provided to support this subdiffusional 
behavior. It would be particularly interesting to try to examine such issues from a rigorous point of 
view, also considering how the relevant subdiffusive scalings may depend on the form (and distribution) 
of the random perturbations. Would similar phenomena also arise in the presence of nonlinear random 
perturbations? Could the integrability of the underlying nonlinear model crucially alter the dynamics 
(e.g., what would happen to the above results if the DNLS was changed to the AL-NLS model)? These 
are only some of the numerous questions that this novel direction raises. 



4. Klein Gordon Models in Mechanical and Electrical Systems 

4. 1 Recent Developments 

Interestingly, in recent years, there has been a parallel stream of developments examining the dynamics 
of the localized modes (discrete breathers — DBs — ) which are in some sense the analogs of the stand- 
ing waves of the DNLS equation considered in the previous section. Some of the fundamental tools 
that enabled the consideration of the existenc e and stability th eory of the discrete breathers and multi- 
breathers were seeded in the original w ork of lAhn et al.\ (1200 lb and subsequently were generalized e.g. 
in that of lMacK av & Sepulchre (2002). Similar results about the stability of such states had been ob- 
tained earlier in lArchilla et al.\ (120031) . using the notion of the so-called Aubry band theory developed in 
Aubrvl 01997b , however the results of the two methods can be shown to be equivalent. In the presenta- 
tion that follo ws, we will be using the formulation o f MacKay and collaborators, as was adapted to the 
KG setting by iKoukoulovannis & Kevrekidis d2009l) . This has the benefit that not only does it allow to 
quantify the existence conditions of the relevant structures, but it also offers tantalizing connections of 
the KG setting with the previously explored DNLS one. It should also be noted that in this setting, as 
well, the relevant results have not been restrained to one-dimensional or just square lattices, but have, 
in fact, been generalized to non-squar e lattices, such as hexagonal and hon eyc omb ones in the works 



of IKoukoulovannis & MacKayl (120051) . IKoukoulovannis & Kourakisl (12.009ft and IKoukoulovannis et al. 
d2009h . 



A short outline of the relevant methodology is as follows. Assume that the system Hamiltonian is of 
the form 



H=H + EH l = £ -pj + V( Xi ) )+- £ (x i+1 - Xi ) 



(4.1) 



where the part of Hq is associated with an unperturbed single oscillator, and H\ denotes the perturbing 
part due to the coupling with the neighboring oscillators. Considering a number of central excited 
oscillators, one can write their solution x(t) into the action-angle form x{t) = £" =0 A, ! (7)cos(nw), where 
J, w are the action-angle variables. Then, over the n excited oscillators, one can define the effective 
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Hamiltonian 



H ef£ = Ho(Ii) + e(#i) i = 1 ...n, 



over the canonical variables: 



# = wo 

01 = w\ — Wq 

(j>2 = W2 — W\ 



si = J + ...+J„ 

h =Jl+...+J n 
h=h+...+J„ 



(4.2) 



(4.3) 



W„-W„-l 



In — J n 



and with 



(H 1 ) = ^H 1 dt, 

where the integration is along the unperturbed periodic orbit. 

One of the key results of the theory of MacKay and collaborators is that the critical points of 
this effective Hamiltonian yield the DB solutions of the model. In the KG setting discussed above, 
Koukouloyannis & Kevrekidis d2009l) derived the effective Hamiltonian in the form: 

1 

2 



= l!=l 



(4.4) 



and hence obtained the persistence conditions through = 0. Remarkably, in this setting, similarly 
to DNLS, it can again be proved that the only available multi-pulse (i.e., multi-breather) solutions are 
the ones with relative phase among the excited sites of or n. 

Subsequently, one can examine the stability of these DB solutions, which is controlled by the Hes- 
sian of the effective Hamiltonian presented above. In particular, E = JD 2 H eS where the symplectic 



matrix reads J 
yields: 



E 



O 
I 



I 

O 



Once again, an explicit calculation of the relevant stability matrix 



A 


B 




D 







Co + eCi 





-e 




£ d 2 (m) 






d<j>id(j)j 


d 2 H 






dlilj 


+ dlidlj 


d(j)jdlj , 



(4.5) 



E 



For the special case of relative phases 0, = 0, or </>, = %, the matrix simplifies considerably acquiring the 
form: 

O B \ _ / O eBj 
C O ) ~ \ Co + eCi O 

which, subsequently, if we consider only the dominant eigenvalue contributions (where X 2 is of < 
yields the relevant squared eigenvalue matrix in the form: 

/ V\ -fx 
-h 2/ 2 -h 



(4.6) 



d(0 d(0 

B 1 .Co--^ 7 Z--- 57 



-fn-\ 






2/„-i 

-fn 



~fn-\ 
Vn ) 



(4.7) 




FIG. 2. (Color online) A sampler of breather solutions of the discrete no nlinear Klein-Gordon equation in one and two spatial 
lattice dimensions with a (soft) potential of the type used in the work of Koukouloyannis & Kevrekidis 12009). The left panel 
shows an out-of-phase two-site mode, while the middle an out-of-phase three-site mode. The right panel shows a two-dimensional 
hexagonal lattice generalization of these mode s in the form of a discrete vort e x brea ther (2jt/3 phase difference between the three 
adjacent excited oscillators), per the work of Koukouloyannis & Kourakis 12009). All three structures are stable close to the 
anti-continuum limit. 



In this expression (0 = dHo/dJ denotes the frequency, while /,■ = f((j>i) = (l/2)^~ = [ n 2 A^cos(«^,). 
This leads to the characteristic exponents (i.e., effective eigenvalues) of the DB in the form: 

(j ±; = ±V^I7+^(e 3/2 ) i = l...n, (4.8) 

where Xn = ~^j z i an d Zj are the eigenvalues of the matrix Z defined above. Remarkably, the count 
of unstable eigenvalues comes out to be identical to the DNLS case, for the focusing or soft nonlin- 
earities which satisfy e^j < 0. That is, the number of unstable eigendirections will be precisely equal 
to the number of nearest neighbors which are in phase between them. Lastly, if the sign of the above 
product e changes, then the conclusion indicates that the number of unstable eigendirections is pre- 
cisely equal to the number of out-of-phase nearest neighbors. Some typical examples of the relevant 
configurations in one- and two-dimensional settings are shown in Fig. [2] 



4.2 Future Perspectives 



While the general Hamiltonian framework of discrete breathers in Klein-Gordon lattices is gradually 
being explored at a fundamental level regarding the existence and stability properties that it entails for 
its coherent structures, numerous new applications are emerging in this field which require considerable 
additional modeling efforts and which bear sig nificant differences fr om this general framework. 

One such example is the recent work of ICuevas et al. where the dynamics of discrete 

breathers were examined in arrays of coupled pendula, in the more realistic experimental setting of forc- 
ing and damping. In the Hamiltonian setting the relevant model for this system is simply the discrete 
sine-Gordon equation, where the linear coupling stems from the torsional springs, and the sinusoidal 
nonlinearity from the gravitationa lly induced torsion. H owever, the actual forced-damped array is more 
accurately described [as argued in lCuevas et al.\ (12009b 1 by: 



2 P y ji ? 

ii„ + co sin(t/„) - jA 2 u„ = -yii„ + y4 2 M« - Fco d cos(co d t)cos(u n ) 



(4.9) 



where, in addition to the intrinsic frequency (Oq of each pendulum and their torsional coupling of strength 
j3, terms for the damping of each pendulum (characterized by y), the viscous damping due to their 
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coupling (characterized by 72) an d the strength F of their parametric drive (characterized also by the 
frequency (O c j) have been included; / denotes the pendulum's mo ment of inertia. In this context, a 
number of interesting features have been highlighted in the work of ICuevas et al. d2009l) . including the 
apparent exchange of stability between the onsite (i.e., site-centered) breather, and the inter-site (i.e., 
bond-centered) form. Both of these waveforms were experimentally observed, each in their range of 
stability, illustrating that there is a separatrix between the respective ranges. The existence of such a 
separatrix connects partially also to the discussion of the previous section about the Peierls-Nabarro 
barriers, and their disappearance, as well as of the latter leadi ng to mobility of the localized modes. The 
mobility of discrete breathers was, in fact, also observed in ICuevas et al.\ (120091) both experimentally 
and numerically. 

Nevertheless, there are still numerous topics within this theme that remain partially unexplored in 
this more elaborate, yet more realistic class of models. On the one hand, it would be useful to explore 
some of the instabilities giving rise to coherent structures therein, such as the modulation al instability. 
Experiments already show the presence of such phenomena as analyzed by English] (120091) . but it would 
be useful if a quantitative description thereof could be provided. Another interesting direction worth 
expanding towards is that of multi-pulses. From the general theory of the previous subsection, an 
understanding is emerging about which multi-pulse solutions (and under what phase differences) should 
be expected to be stable/dynamically robust. Developing a similar framework of understanding for the 
forced-damped setting initially from a computational perspective (and possibly also analytically at a 
later stage) would certainly be desirable. 

A similar driven damped se tting has recently been e xplored in the conte xt of driven and damped 
electrical lattices in the work of lSato et al.\ (12007b and of English et al.\ ( 12008b . There the experimental 
system is a bi-inductive discrete electrical transmission line, which is nonlinear due to the nonlinear 
dependence of the capacitance of the diode on the voltage and which leads to the following dynamical 
model: 



G>qU„ ■ 



01 



22o 3(2 



CO, Aiu n 



co? 2 cor 1 dV d 

~A 2 uZ-\ KA 2 u„-\ f- 

2Qo " 3Ql " R dt 



1 



1 



(4.10) 



Various electrical constants are denoted by their respective symbols (R for resistance, Co for capacitance 
and Qq for charge), while the external driver is given by V e /. Some of the very interesting features 
observed in this system include how its discrete breathers could interact with impurities (being e.g. 
seeded by or attracted to a capacitor, while they could be destroyed or repelled by inductors and/or 
resistors). Also, it was found that depending on the driver amplitude and especially the driver frequency, 
the system would spontaneously lock into 1- or 2- or multi-breather states within a finite size (e.g. with 
24 nodes) experimental lattice. It would be extremely interesting to try to explore the existence and 
the stability of localized modes within such transmission line systems and to explore the origin of 
the observed pattern selection and dynamics. To the best of our knowledge, very little has been done 
theoretically in this direction for such more realistic systems. 

Interestingly, electrical lattices constitute, arguably, the simplest realization of one of the focal points 
of interest in the last few years (which is lik ely to be one of the major areas of development for the next 
decade), namely of nonlinear metamaterials lKozvrev & van der Weidd 020081) . Such media have simul- 
taneously negative electric permittivity and magnetic permeability (so that their root square product, i.e., 
the speed of light propagation is real); however, because of these negative values, they feature a num- 



Nonlinear Waves in Lattices: Past, Present, Future 



21 of|35] 



ber of unusual effects, including negative refraction, opposite signs of group and phase velocity (i.e., 
wavefronts and wavepackets move in opposite directions), inverse Doppler effec t, backward Cerenkov 
radiation, etc. A nice introduction to this field can be found e.g. in the work o f Marques et al. d2008h 
[see also for a review with a more nonlinear slant the presentation of Lichtnitser & ShalaevI ( 2008l) l. 
One of the simplest (if not the simplest) lattice that produces such backward waves can be generated 
by feeding a ladde r network with alterna ting shunt-connected inductors and series capacitors [as shown 
e.g. in Fig. 3.1 of iMarques et al.\ d2008l) l. However, recent work in the realm of nonlinear waves has 
mostly focused on the slightly more elaborate setting of magnetic metamaterials and, in particular, 
lattices of split-ring resonators. The resonators are electrically equivalent to an RLC oscillator, with 
self-inductance L, ohmic resistance R and capacitance C and become nonlinear because of the gaps that 
they possess which are filled with a nonlinear dielectric medium (whose permittivity depends on the 
intensity of the electric field E). The relevant mathematical description of the dynamics of the charge 
stored on the capacitor of the «-th split-ring resonator (SRR) of such a SRR lattice, coupled to the current 
circulating through it can be obtained as follows: 



dQ„ 
dt 

Ldl„ 
dt 



= In 



-IR„+M 



dl n 



+ 1 



dt 



din-] 

dt 



-/(&.)+*. 



(4.11) 
(4.12) 



Here, <§ is the electromotive force (i.e., external drive) induced in each SRR due to the applied field, 
while f(Q n ) bears the nonlinearity of the model since the field-dependent permittivity gives rise to a field 
dependent capacitance, which in turn leads to its nonlinear dependence on the charge. / can be reason- 



ably a pproximated by an odd power series in Q„, as is illustrated e.g. in the work of [Lazarides et al. 



This formulation of, o nce again, a damped and driven electrical lattice has been used in the 
work of lLazarides et al.\ (120061) to predict the existence of dissipative breathers (as well as Hamiltonian 
ones, in the absence of the driving and damping) under an AC (alternating current) electromotive force. 
This work initiated a path towards a more sy stematic exploration of such structures, which was con- 
tinued in the work of Elefthe riou et al.\ d2009l) . examining such structures in both 1- and 2-dimensional 



sett ings, that of iMaluckov et al. I d2008h for different types of nonlinearities, of Lazarides et al. 



and Eleftheriou et al. (120091) for surface variants of these discrete breathers and of Molina et al.\ 



(2008) 



2009) 



for binary metamaterials. More recently, it should be mentione d that further "frequ ency gap" regions 
in such metamaterials have been explored e.g. in the work of iTsitsas et al. (1201 Oh . where the signs 
of the magnetic permeability and dielectric permittivity are opposite; this regime has been predicted 
to al so possess breathers due to its approximate description by the so-called short-pulse equation [see 



e.g., ISchaeffer & Wayne! (120041) 1. However, it should be mentioned that in most of these settings, the 



relevant models are derived, and subsequently are analyzed in a purely numerical form. Both a theoret- 
ical/mathematical analysis of the types of structures that are possible, and experimental realizations of 
the relevant settings are presently lagging behind and would be a natural direction for future work. 



5. FPU-type Models in Granular Crystals and Dusty Plasmas 

5.1 Recent Developments & Future Perspectives 

FPU models have a time-honored h istory of more than ha lf a century, one that has been captured in 
bo th special volum es such as that of ICampbell et al.\ (120051) . as well as in recent books such as the one 
bv lGalavottil (12008b . Despite the interesting recent developments in this area, including the ^-breathers 
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of iFlach et al\ ( 12005b (which may be closely related to the original FPU paradox), here we will not 
focus on this intensely-studied theme. Instead, we will present some particular recent directions within 
the same general class of models which feature predominantly intersite nonlinear interactions within a 
lattice setting. 

The principal example of the above class that we will consider is that of the so-called "granular 
crystals", namely arrangements of densely packed spheres, which upon suitable excitations interact 
elastically through the so-called Hertzian contacts, and result in the propagation of compression waves 
(or, as we will see, the localization of suitable breather-like states) in these media. These chains have 
drawn considerable attention during the past few years. This broad interest has emerged, to a consid- 
erable extent, due to the wealth of available material types/sizes and the ability to tune their dynamic 
respon se to encompass lin e ar, weakly nonlin e ar, and strongly non linea r regimes as s hown in the works 
of e.g. iNesterenkof d200ll) . ISen et al\ d2008l) . IDaraio et al\ 00061) and ICoste et al\ dl997h . Such flex- 



ibility makes them perfect candidates for m any engineering app l ications, incl u ding shock and energy 
abs orbing layers p r opose d e.g. in the works o flDaraio et al\ (120061). Hong ] d2005l) . lFraternali et al\ (1201 On 



and Donev & Sen (120061) : a ctuating devices [ s ee e.g . iKhatri et al. and sound scramblers 

e.g. IDaraio et al.\ (120051) and lNesterenko et al.\ (120051) 1. 

As an aside, we should note here that in an area such as that of granular crystals, which from the 
perspective of mathematical results is still at an early stage (most of the focus has been on numerical 
investigations and experimental results), it is not easy to distinguish between recent developments and 
future perspectives. For this reason, (and differently than in the previous sections) we will present both 
of these in a unified format in what follows. 

The model used to represent the simplest form of such granular systems consists of a one di- 
mensional (ID ) chain of spherical beads regulated by Hertzian contact interaction potentials, see e.g. 
Johnson! (120041) . For this class of models, a nonlinear wave theory, supp orting nearly - comp act solitary 
waves was derived for all structured homogene ous [as discussed e .g. in iNesterenkol (120011) 1 and later 
also for heterogeneous periodic [as illustrated in lPorter et al.\ (120081) 1 materials showing a highly nonlin- 
ear force (/^-displacement (8) response dictated by the intrinsically nonlinear potential of interaction 
between its fundamental components. These general nonlinear spring-type lattice interactions read: 



-V'(u r 



)=A8 P , 8 = max(0,w,-_i — it,-), 



-V'iui-ll^+V'illi+y-lli) (5.1) 



where A depends on the material properties and p is the nonlinear exponent of the contact interaction 
(with p > 1); p = 3/2 for spheres (Hertzian interactions). By granular matter, we mean an aggregate 
of particles in elastic contact with ea c h other, preferably in li near or network shaped arrangem ents 



as discussed e.g. i n iGoddardl dl990» . iGilles & Costel d200lh . iGoldenberg & Goldhirschl d2005l) and 



Hostler & Brennenl (120051) - In addition to the nonlinear contact interaction, and related purely to the 
particle geometry, another unusual feature of the granular state is the zero tensile strength, which in- 
troduces an additional nonlinearity (asymmetric potential) to the overall response. In the absence of 
static precompression, these properties result in an essentially nonlinear force which leads to interesting 
implications on the d ynamics and a vanishing sound speed co = (so-called "sonic vacuum"); see e.g. 
Vakakis et al\ (|2008). This makes the linear and weakly nonlinear continuum approaches based on the 
Korteweg-de Vries (KdV) equation invalid and places granular crystals in a strongly nonlinear regime 
of wave dynamics. This feature uncovers a significant potential for relevant applications. In particular, it 
supports a new type of nearly compact highly tunable solitary waves that have been experi mentally and 
numer icall y observed for ID Hertzian granular systems as summarized e.g. in the reviews of lNesterenkol 
(I200lb and lSen et al.\ ( 12008b . While spherical bead crystals have b een mostly studied, s imilar responses 
can be obtained from cylindrical and elliptical particles [see e.g. IDaraio et al.\ d2009l) l. fibrous layers 
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[see e.g. Lambert & Tesand 19841) 1 and foams [see e.g. lSen et al.\ ( 120081) 1. 

The prototypical setup of Hertzian contacts describing elastic interactions among grains of a unary 
(i.e., monoatomic) lattice in the absence of any precompression (i.e., in the so-called sonic-vacuum 
regime) can be written as: 



u n = [m„_i — u„\ + — [u n — U n+ \\ + . 



(5.2) 



Here, the [•] + notation indicates that that the quantity in the bracket is only evaluated if positive, while 
it is set to 0, if negative. Reformulating the problem based on the strain variant of the equation for 

r n = u„- \ — u n , we obtain 



■2[rn]t 



(5.3) 



If we now seek traveling wave solutions, as suggested by the above mentioned numerical and experi- 
mental work (as well as the nea r-co ntinuum asymptot ics used) by numerous authors [see e.g. the pre- 
sentations of Nesterenkol ( 2001 ) and Sen et al. (2008)], then r n = (j)(x) = (j)(n — ct). This leads to an 
advance-delay equation the solution of which, from a rigorous viewpo int, "remains an open question 
and apparently a delicate one" as indicated by English & Pegol (120051) . This is because, although the 
latter problem can be reformulated as an iterative scheme of the form: 



<Pm 0) 



(5.4) 



with A(x) = [1 — |x|]+, the proof of existence of a localized solution thereof for any p > 1 is presently 
missing. Nevertheless, from a numerical pe rspective, the scheme o f E q. (15 .41) has been shown to con- 
verge to the desired solutions in the works of English & Pegol ( 2005 ) an d lAhnert & Pikovskvl(l2009l) : for 
an example of the relevant traveling wave in the Hertzian case, see Fig. [3] It should also be mentioned 
that should such a proof be available, it would immediately imply that 



r (x + 1 ) < sup vG { _ lA] rP(x+\-y) J_A{x)=r p {x), 



(5.5) 



which, in turn, suggests the double exponential law r(x + n) ^ r(x) p " for the decay of the solution tail. 
It is this extremely fast decay which leads the continuum analog of the solution to sustain compactly 
supported structures. However, it should be noted that this latter statement has not been rigorously 
proven either (to the best of our knowledge). In fact, there are multiple continuum models that have 
been proposed for the system dynamics. The most standard of them stems from the long wavelength 
approx imation (LWA) applied to the displacement pr oblem bv iNesterenkol (feOOll) [see also Coste et al. 



1997 )1. while a more recent approach in the work of Ahnert & PikovskyT i 20091) suggests that the LWA 



be applied to the equation for the strains. The former approach leads to the partial differential equation 
(PDE), for the strain variables, of the form: 



while the latter establishes: 



(r p ) xx + —[ (r p ) xxxx + 



n{n — 1) 



- 2 rl) 



(5.6) 



r n = (r p ) xx + —(r p ) xxxx . 



(5.7) 
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FIG. 3. (Color online) A typical example of the spatial profile of the positions (blue solid line connecting the lattice nodes) and 
of the velocities (red dashed line connecting the lattice nodes) of an exact (nearly-compact, see inset) traveling wave solution of 
the monomer granular chain, obtained through the solution of Eq. 45.3t . The space (n) - time (t) evolution of this wave with its 
prescribed speed of unity is shown through the contour map of the strain r n (t). 



In fact, from a mathematical perspective, one can go even further and suggest that since the above models 
are singularly perturbed (and may suffer relevant pathologies, which may also affect their computational 
implementation), it may be relevant to s uggest the consi derat ion of regularize d variants of these PDEs, 
such as the one proposed by the work of lRosenaul (119861) and lRosenaul (120031) as: 



r t t = {r p ) xx + —r xxtt . 



(5.8) 



A systematic examination of the relevant models and a rigorous proof of the decay properties of its as- 
sociated solutions (and how they connect to those of the discrete problem from which the approximation 
stems) would be a particularly interesting theme for future studies. 

A recent direction of intense interest within this theme of granul ar crystals has conc erne d the exami 



nation of progressively more heterogeneous settings. In the work of lPorter et al.\ d2008l) and lPorter et al. 



d2009h . dimers of different types were considered (e.g., 1 steel bead-1 teflon bead, or 2 steel-1 teflon, 
or 5 steel-1 teflon etc., as well as similarly with other materials such as brass, rubber, glass and 
nylon), as well as of trimers (1 steel bead, 1 brass bead, 1 teflon bead or 1 steel, 1 glass, 1 nylon, etc.). 
From a theoretical/mathematical perspective, one of the fundamental contributions of this work for the 
consideration of 



miiij = (wj - ujf -{uj- Wj_]) k , 
m 2 Wj = (u j+l -Wj) k -(wj-Uj) k , 



(5.9) 
(5.10) 



(i.e., the reduced dimer granular setting) was the proposition of a modified Taylor expansion: 

w = X (u + b x Du x + b 2 D 2 u xx + b 3 D 3 u xxx + b 4 D 4 u 4x + ...). (5.11) 

This enabled the effective "homogenization" of the ordinary differential equations into a "dimer contin- 
uum medium", provided that the relevant consistency [between Eqs. d5.91 l and (15.1 OH conditions A = 1, 
b\ = \,b% = mi /(mi +"12), bj, = (2mi — m^) / (3 (mi -\-rn2)), b 4 = mi (m\ — 11111112 (2>{nii +1112)^), 

etc. are satisfied. Notice that for mi = m2, these yield the regular Taylor expansion, as expected. This 
long-wavelength approach for the dimer problem produced a PDE of the type of Eq. (15.61 , whose so- 
lutions were comp ared favorably to nume rical computations and experimental results in the works of 
Porter et al.\ d2008l) and lPorter etal\ d2009l) . 
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On the other hand, the consideration of dimers enabled for the first time another possibility that has 
been very recently explored more systematically. In particular, in the presence of precompression, the 
dimer model becomes: 



litti =A[8 Q + u ( -i - uj} 3/2 - A[<5 + ui - u i+ i] 



3/2 



(5.12) 



where m.2i+\ = M and m^x = m for all i £ Z (for same material, different radii Rj beads in the simplest 
realization). Then, at the end of the first Brillouin zone, i.e., at k = where a = Rj +Rj+i — So 
is the equilibrium distance between two adjacent beads, the linear spectrum possesses a gap between 
the upper cutoff of the acoustic branch, (Q\ 

-J 2K2 1 m, where K2 

an approximate K% — K$ — K4 chain of the form: 



\J2K2 /M, and the lower cutoff of the optical one, (xh 

\I2 

fAaft'". This gap and the corresponding reformulation of the granular chain as 



niM = K 2 (u i+ \ - 2iij + Ui-i) +K3 ((«i+i - ui) 2 - (m,_i - m,-) 2 ) 
+ K 4 ((m /+1 - ui) 3 + (m,-_i - ui) 3 ) 



(5.13) 



with K3 = — |A5q ^ 2 and K4 = ^A5 J/z enables the adaptation of the findings of lHuang & Hul d 1 998b 



-3/2 



demonstrating that since > |, discrete gap breather modes will bifurcate from the optical band, in 
fact, due to the modulational instability of the band edge. The relevant nonlinear states were recently 
experimentally and numerically established in the work of iBoechler et al. d2009l). Furthermore, s uch 
breathing excitations were also established in the presence of impurities by lTheocharis et al.\ (120091) . as 
a nonlinear continuation of the linear modes induced by the presence of such defects. In that setting 
the relevant nonlinear impurity modes were obtained both in the setting of one defect, as well as in that 
of a defect "double well" (with two defects separated by one lattice node). Similar excitations but of a 
transient form, due to the absence of precompression (and hence of an underlyin g linear regime), were 
also recently observed experimentally (and discussed theoretically) in the work of ljob et al.\ (120091) . 

It is interesting to also note that once the role of individual defects is starting to be explored, it is 
natural to examine the progressive inclusion of disorder within suc h granular chains. A few steps have 



been already taken in that direction. In particular, in the work of iPonson et al\ (120091) . a systematic 



approach towards increasing randomness within the lattice has been used, by examining a dimer setting, 
and progressively "flipping" the order of the masses in the dimer in order to create defects. In that case, 
a transition between the propagation of waves in low-disorder regimes and the degradation of the wave 
in high-disorder lattices was clearly illustrated. On the other hand, the incorporation of disorder can be 
used to achieve a degr ee of optimality in granular networks. An example in this direction is the work of 
Fraternali et al. where upon a given selection of masses and elastic properties thereof, a genetic 

algorithm approach was used towards the optimization (minimization) of the propagated force at the 
end of the chain. This was a first step towards the optimal construction of "granular protectors". 

The above investigations are indicative of the intense recent activity in the field of granular crys- 
tals. Yet, there is a wide range of topics that require additional theoretical, numerical and experimental 
investigation. In the one-dimensional setting, in addition to exploring the existence (and stability) of 
traveling waves from a rigorous viewpoint not only in the monomer, but also in the "multimer" (i.e., 
dimer, trimer etc.) setting, it would also be of interest to explore more systematically the role of disor- 
der. In particular, the setting of granular chains might be an ideal testbed where the recent developments 
of subdiffusive wavepacket spreading and possible "destruction" of Anderson localization in nonlinear 
lattices could be explored controllably. This is because this setting can be made, on demand, nearly 
linear, weakly or strongly (or even completely) nonlinear. Furthermore, while the above ideas discussed 
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earlier bas ed on the works of | Kopidakis et al.\ d2008b . iFlach et al. d2009l) . iPikovskv & Shepelvanskvl 



(2008) and Skokos et al. ( 20091) . were established in the context of nonlinear models with onsite po 
tentials such as the DNLS, it would be interesting to explore them in models within the FPU class i.e., 
with inter-site interactions, and especially ones such as the granular chain which could in principle be 
realized experimentally in a controllable fashion. On the other hand, another important observation 
in many of the above men tioned experimen tal works in granular chains [see also the relevant reviews 
of iNesterenkol (1200 lb and ISen et al\ d2008l) l is the key role of dissipative perturba tions in decreasing 
the amplitude of the propagating [or localized breathing as in lBoechler et al.\ d2009b l waves. A widely 
applicable, quantitatively accurate modeling from first principles of such mechanism s is still miss i ng de- 
spite r elevant recent efforts an d interesting corresponding results e.g. in the works of lRosas et al.\ (120071) 
and of Carretero et al. d2009h . Another direction which still remains relatively pristine is that of higher 
dimensional structures. While some of the linear or near-linear behavior has been explored therein in 
both ordered and disordered granular chains, much less is known about the structure and propagation 
(as well as the localization) of nonlinear waves within such higher-dimensional granular crystals. This 
is a theme which, considering the activity and progress recently spurred in the one-dimensional setting, 
is likely to emerge as one of the most lively areas of investigation within the next decade. 

Finally, we wanted to add a brief comment on FPU-type lattice applications in dusty plasmas and 
the potential emergence of nonlinear excitations in these settings. In rece nt works such as those of 
Koukotilovannis & Ko urakis (2007), Koukoulovannis & Kou rakis (2009) and Kourakis & Shukla ( 20061) . 



the setting of dusty plasmas has been proposed as a canonical realization of nonlinear Klein-Gordon lat- 
tices. The coupling between adjacent nodes/dust grains has been proposed to be linear, while the sheath 
environment has been argued to create a nonlinear onsite potential landscape that drives the relevant 
transverse dynamics considered in these works. In this context, both one- and two-dimensional discrete 
breathers of the type explored in the previous section for Klein-Gordon lattices have been identified 
and their stability has been accordingly examined. However, here, we would like to foc us the atten- 
tion instead on the broader context of e.g. the presentation of iKourakis & Shuklal d2006b . There both 
longitudinal and transverse vibrations of the relevant dust p articles were exam ined. Although the latter 
are, in principle, experimentally possible as in the work of IZafiu et al.\ (1200 lb . there is also a particu- 
lar interest in the former, as they incorporate the nonlinear potential of interaction across sites which 
can be accurately approximated by an exponentially dependent on the inter-particle separation Yukawa- 
type landscape. Both the linear and the n onlinear implication s of such lattices have been explored in 



some detail as e.g. in the recent review of iMelzer et a l. (2008). However, the emergence of nonlinear 



waveforms and especially the existence, stability and mobility of discrete breathers in these intersite but 
non-nearest-neighbor (in fact, the coupling is "all-to-all") lattices have received limited attention. We 
feel that this field is also presently at a ripe phase, where the use of the tools and technology of local- 
ized modes and their stability and dynamics could have a substantial impact in experimental, as well as 
numerical investigations. 



6. Instead of an Epilogue 

The field of nonlinear dynamical lattices has seen a tremendous explosion of interest associated with it 
over the past 20 years. Although the first anharmonic lattice investigations were initiated more than 50 
years ago by Fermi, Pasta and Ulam, and part of these led to the remarkable theory of solitary waves 
and solitons which has been pervasive of various areas of science for more than the past 40 years, the 
appreciation of the idiosyncrasies of the lattices in comparison to their continuum siblings is far more 
recent. The seminal studies on intrinsically localized modes led to remarkable breakthroughs that fu- 
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eled research in areas as diverse as nonlinear optics and optical waveguides, Bose-Einstein condensates 
in optical lattices, micromechanical systems, pendulum arrays, electrical lattices, DNA denaturation, 
Josephson junction arrays, split ring resonators and metamaterials, granular crystals, dusty plasmas, 
metallic uranium and ion i c sodi um iodide, and possibly even in sheets of graphene [see e.g., the recent 



work of lYamayose et al.\ (120071) 1. Alongside the development of the areas of application, the particu- 
larities of each of these settings have forced not only a development of the general mathematical ideas 
and of fundamental conditions under which the relevant excitations should be expected to exist and be 
robustly preserved by the system dynamics, but have also inspired the more specific adaptation of these 
ideas and of the relevant modeling to the special settings of each particular application. They have also 
motivated the emergence of a diverse array of relevant excitations, ranging from solitons to discrete 
breathers, to multi-peaked structures or even to compactly (or nearly compactly) supported ones. It is 
clear that after the initial decade of enthusiasm and exploration, and the second one of broad and far 
reaching applicability, the field is now entering a new decade of maturity, but where significant chal- 
lenges lie ahead. These concern the application of the relevant ideas to very recently emerging fields 
(materials and metamaterials, granular media, plasmas), their exploration in more realistic and espe- 
cially higher dimensional settings and the establishment of more firm connections/interplays with other 
areas of physics, such as non-equilibrium thermodynamics, or the study of disordered media and of 
wave phenomena within them. What the resulting 20/20 vision of nonlinear waves in lattices will be 
remains to be seen. . . 
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